Dark solitons and vortices in PT-symmetric nonlinear media: 
from spontaneous symmetry breaking to nonlinear VT phase transitions 
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We consider the nonlinear analogues of Parity-Time (VT) symmetric linear systems exhibiting defocusing 
nonlinearities. We study the ground state and excited states (dark solitons and vortices) of the system and report 
the following remarkable features. For relatively weak values of the parameter e controlling the strength of the 
^T-symmetric potential, excited states undergo (analytically tractable) spontaneous symmetry breaking; as e 
is further increased, the ground state and first excited state, as well as branches of higher multi-soliton (multi- 
vortex) states, collide in pairs and disappear in blue-sky bifurcations, in a way which is strongly reminiscent 
of the linear PT-phase transition — thus termed the nonlinear PT-phase transition. Past this critical point, 
initialization of, e.g., the former ground state leads to spontaneously emerging "soliton (vortex) sprinklers". 



Introduction. Over the past decade, and since its original 
inception 0112], the theme of VT- symmetric Hamiltonians 
has gained considerable momentum in the physics and ap- 
plied mathematics communities. Such systems, respecting 
both Parity (V) and Time-reversal (T) symmetries — still ex- 
hibiting real spectra while non-Hermitian — provided an in- 
triguing alternative to standard Hermitian quantum mechan- 
ics. Note that for a standard Schrodinger type Hamiltonian 
with a generally complex potential U, the VT symmetry dic- 
tates that the potential satisfies the condition U (x) = U*(—x) 
[where (•)* stands for complex conjugation]. 

Despite the theoretical appeal of such models, it was only 
recently shown [3] that optics could be an ideal playground 
for the physical/experimental realization of systems featur- 
ing the VT symmetry. However, this also added another ele- 
ment in the interplay, namely nonlinearity. In that context, the 
considerations of Ref . yfl extended from bright and gap soli- 
tons to linear (Floquet-Bloch) eigenmodes in periodic poten- 
tials, examining how these coherent structures are affected by 
the genuinely complex, yet T^T-symmetric potentials. More 
recently, experimental results were reported both in nonlin- 
ear optical systems flBl] and electronic analogs thereof ((J. 
These, in turn, have triggered a wide range of theoretical stud- 
ies on nonlinear lattices with either linear fF7l— Tl~3ll or nonlin- 
ear lfl6l - [T8h PT-symmetric potentials and, more recently, on 
harmonic PT-symmetric potentials 11911 . 

While the above volume of work has examined numerous 
features extending from bright solitons to defect modes, and 
from gap solitons to PT-lattices, the consideration of defo- 
cusing nonlinearities, and especially of dark solitons has been 
extremely limited (see, e.g., Refs. Il20ll ). Little attention (and 
again chiefly in the focusing nonlinearity case JH]) has also 
been paid to PT-symmetric systems in higher dimensional 
settings and the corresponding interplay with nonlinear states 
such as vortices. In the present work, we study systems with 
PT-symmetric Hamiltonians exhibiting defocusing nonlin- 
earities, and focus on the existence, stability and dynamical 
properties of the ground state and excited states, i.e., dark soli- 
tons and vortices. Our main findings for a prototypical VT- 
symmetric potential, which is harmonic in its real part and 



has a localized imaginary part (parametrized by an amplitude 
parameter e) are summarized as follows: 1) dark solitons are 
shown to be subject to spontaneous symmetry-breaking (SSB) 
instabilities for small e\ 2) for higher values of e, the ground 
state and the first excited state (single dark soliton), as well 
as pairwise — e.g., 2nd and 3rd, 4th and 5th etc. — higher ex- 
cited states (respective multiple dark soliton states) are subject 
to a nonlinear analogue of the PT-phase transition, collid- 
ing and disappearing in a set of blue- sky bifurcations; 3) be- 
yond this critical point, the system acts as a soliton sprinkler, 
spontaneously emitting dark multi-soliton structures. 4) All 
of these features have direct counterparts for vortices in two- 
dimensional settings, illustrating the generic nature of these 
findings. 

Fundamental States: Ground State and Dark Soliton. Our 
model, motivated by the above nonlinear optical considera- 
tions (but also by ones pertinent to nonlinear phenomena in 
Bose-Einstein condensates (BECs) i2lll ). will be, for the one- 
dimensional (ID) setting, as follows: 

id t u = ~\d 2 x u + \u\ 2 u + [V(x) + iW{x)]u, (1) 

where u is the complex electric field envelope (or the macro- 
scopic wavefunction in BECs), t denotes the propagation dis- 
tance (or time in BECs) and x is the transverse direction. For 
a PT-symmetric Hamiltonian, the real and imaginary parts 
of the potential must satisfy V(x) = V(—x) and W{x) = 
—W(—x). Below we focus on the case of a real parabolic 
potential, V(x) = (l/2)Q 2 x 2 , modeling the transverse distri- 
bution of the refractive index (or the external trap in BECs), 
while the imaginary part W(x) is considered to be an odd, 
localized function of spatial width <C Q -1 . A generalization 
of this model will be studied below in two-dimensions (2D), 
with V = (l/2)n 2 (x 2 + y 2 ) and W = e(x + yJe-^+f 2 )/ 4 . 

We now analyze the fundamental states (namely the ground 
state and the first excited — single dark soliton — state of the 
system) shown in Fig. [T] We seek stationary solutions of 
Eq. (Q]) in the form u = Ub(x) exp(— ifit), where fi is the prop- 
agation constant (or the chemical potential in BECs). For a 
sufficiently small imaginary potential, W(x) = eW(x) [with 



2 



max{|VK(x)|} = 0(1)], where £«1, and when the inverse 
width ft -1 of V(x) is sufficiently large, ft ~ £, we may find 
— in the Thomas-Fermi (TF) limit — an approximate solution 
of Eq. (Q]). This is of the form = [^f/I + / (x)] exp[i(/>(x)] 9 
where the amplitude and phase f(x) and 4>(x) (considered to 
be small, of order e 2 and e, respectively) are given by: 



/(*) = 



2v^ 



(V + 2W 2 



(j)(x) 



Wdx. (2) 



where W = / VKdx. Contrary to the conservative case 
(e = 0) i2lll . this TF background is characterized by a den- 
sity dip located at the center (x = 0) and a nontrivial phase 
distribution. Both features are shown in the top panels of 
Fig.[TJ where the analytical result is compared with the numer- 
ical. Importantly, a linear stability — Bogoliubov-de Gennes 
(BdG) — analysis (see, e.g., Ref. l25ll ) shows that the back- 
ground U}) (x) is stable against small perturbations. Here we 
notice that the evolution of the power, N = J \u\ 2 dx, is gov- 
erned by the equation dN/dt = 2 j \ui)\ 2 W(x)dx and, thus, 
for U}) even and W(x) odd, it is conserved. This conservation 
also holds for excited states sharing the parity of 1 2x5 1 2 . 

To describe the dynamics of a dark soliton (DS) v(x,i) 
on top of this TF background, we introduce the ansatz u = 
Ub(x)v(x,t) into Eq. (Q} and, after using the scale transfor- 
mations t — >• fit and x —> ^JJix, we obtain the following per- 
turbed nonlinear Schrodinger (NLS) equation: 

id t v + \d xx v + v(l - \v\ 2 ) = /i- 2 P(^), (3) 

where the perturbation [which is of order 0(s 2 )] is given by 



P(v) = (l-\v\ 2 )v (V + 2W 2 )+v x (~V x -2 



(W-i)W) . 



Applying the perturbation theory for dark solitons 12511 . we 
seek a solution of Eq. ([3]) in the form of the dark soliton of the 
unperturbed system (P(v) = 0): v(x, t) = cos <p(t) tanh £ + 
i sirup (t), where £ = cos ip(t) [x — xo(t)], with cp(t) and 
xo (t) being the slowly-varying phase and center of the soli- 
ton. In the adiabatic approximation dxo / dt = sin cp, while 
the perturbation-induced evolution equation for ip reads: 



dip 
~dt 



\d x V- J sech 4 (0 [tanh(0W 2 + WW] dx, (4) 



where we have assumed almost black solitons with cos <p ~ 1. 
This way, for a given W(x), we can derive an equation of mo- 
tion for the soliton center xq . Hereafter, we consider an imag- 
inary potential of the form W(x) = exexp(— x 2 /2) (other 
choices, e.g., W = e sech 2 (x) tanh(x) have led to similar re- 
sults). To examine the stability of the equilibrium at xq = 0, 
we Taylor expand Eq. ©, obtaining to leading order 



d 2 x 



-uj osc x , 



V2 



(5) 



Equation © implies that if the amplitude e of W(x) is less 
than a critical value e$ = ^5/12Q, the soliton performs 
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FIG. 1: (Color online) Fundamental states. Top: Density (left) 
and phase (right) of the numerically obtained TF background [solid 
(blue) line] compared to the prediction of Eqs. ® [dashed (red) line]; 
the inset shows the characteristic density dip induced by W(x) at the 
origin. Bottom: The same properties are shown now for the (single) 
dark soliton. Parameters values are: ji — 3, Q = 0.1 and e — 0.3. 



oscillations in the complex potential with frequency cc; osc ; 
on the other hand, if e > e^J the soliton becomes unsta- 
ble. The above prediction has been confirmed numerically, 
both by means of direct simulations and employing a BdG 
analysis. The latter reveals that the considered stationary 
dark soliton is characterized by the anomalous mode eigen- 
frequency uo a 10,1251], which is real for e < e^J (in this case, 
u a = ^osc)> and it becomes imaginary for e > e$ , thus 
signaling the onset of the SSB instability of the dark soliton 
(which displaces the soliton from the trap center). 

The dependence of uu^ on the amplitude e of the imaginary 
potential W, as found by the BdG analysis, is illustrated in the 
top panel of Fig. [2 As shown in the inset, u\ initially moves 
towards the spectral plane origin, and past the critical point, 
e^cr (cf. vertical line), uo a exits as an imaginary pair of eigen- 
frequencies manifesting the soliton 's exponential instability. 
As shown in the top panel, for small e, the agreement between 
the analytical prediction of Eq. §5§ [(red) line] and the BdG 
numerical result [(blue) circles] is excellent. 

From symmetry breaking to nonlinear VT transitions. For 
larger values of e, the unstable imaginary eigenvalue, makes a 

maximal excursion along the imaginary line and returns to the 

(2) 

origin at a second critical point, Sc r = 0.62, finally collid- 
ing with it. The branch of single soliton solutions disappears 
past this critical point. To better understand how the branch 
ceases to exist, we first observe (bottom panel of Fig. Q]) that 
the density profile of the soliton becomes increasingly shal- 
lower (i.e., more "grey") as e grows and the second critical 
point is approached. This is due to the development of an 
increasingly strong even imaginary part of the solution. Fur- 
thermore, the stable background (ground state) solution Ub(x) 
[cf. Eqs. © and top panel of Fig.[U develops an odd imagi- 
nary part resembling a (progressively darker) grey soliton. Fi- 
(2) 

nally, at e = Sc r , the profiles of these modes become identical 
and disappear in a blue- sky bifurcation through their collision. 
This is shown in the bottom panel of Fig. [2 where the power 
N is shown as a function of e. The top solid (blue) branch 
shows the stable ground state, u^, which ultimately collides 
with the one soliton (first excited) state at e w 0.62 (for fi = 3 
and ft = 0.1). 

Importantly, we have confirmed that the above description 
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FIG. 2: (Color online) Bifurcations/instabilities of nonlinear states. 
Top panel: the dependence of the soliton's squared eigenfrequency, 
cj^, on the amplitude e of W(x)\ the inset shows a detail near the 
critical value siV (vertical line), where it crosses the origin. The 
(blue) circles and the (red) line depict, respectively, numerical (BdG) 
and analytical [cf. Eq. ([5])] results. Bottom panel: the full bifurcation 
diagram (denoting the above instability by a vertical line and) en- 
compassing the pairwise bifurcations/disappearances of the nonlin- 
ear states. Solid (dashed) lines indicate dynamically stable (unstable) 
branches. Here, \i — 3 and Q = 0.1. 



holds also for higher excited states (multiple dark soliton so- 
lutions), as shown in Fig. [2j each pair of the higher excited 
states (2nd with 3rd, 4th with 5th etc.) also disappears in a 
blue- sky bifurcation. A general remark is that higher excited 
states bifurcate for larger values of s. Remarkably, this can 
be thought of as a nonlinear analogue of the VT transition, in 
analogy with the pairwise collisions in Ref. (]]] (see e.g. Fig. 1 
of that reference) for the linear setting ll26ll . 

A relevant question concerns the dynamics of the nonlinear 
waves when subject to these (SSB and blue-sky) bifurcations. 
To answer this, we numerically integrated Eq. (Q]) as shown in 
Fig. [3] In the top panels, we have illustrated the dynamics of 
the DS upon its destabilization at e = . When the SSB is 
manifested, the soliton is either spontaneously ejected towards 
the lossy side (and typically found to localize therein while the 
background grows in amplitude and widens) or moves to the 
gain side executing oscillations thereafter. On the other hand, 

(2) 

past e = Scr , using, as an initial condition the form of the 
TF background (bottom panel of Fig. O, we have found that 
a dark soliton train is spontaneously formed, with an increas- 
ingly larger number of solitons as larger values of e are used 
(i.e., a "soliton sprinkler" emerges). This can be intuitively 
connected to the observation of Fig. [2] that higher excited 
multi-soliton states persist for larger e than lower ones. Again, 
it is typically observed that the solitons are nucleated and stay 
in the vicinity of the global minimum of W(x), which indi- 
cates the "lossy" side of the imaginary potential. 

Two-dimensional Generalizations. The bifurcation of the 
nonlinear structures emerging in 2D follows a similar, but also 
more complex, pattern than in the corresponding ID setting. 
Figure @] depicts the full bifurcation scenario for solutions 
bearing no vortices (the TF background cloud), one to six vor- 
tices, and the dark soliton stripe. As in ID, the TF background 
is stable in all its domain of existence and collides, in a blue- 
sky bifurcation, for a large enough value of e, with an excited 
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FIG. 3: (Color online) Bifurcation-induced dynamics. Top panels: 
manifestations of the SSB destabilization scenarios for an unstable 
dark soliton past e = e^cr . Bottom panels: soliton sprinkler spon- 
taneously leading to two or four solitons from the ground state used 
for e > Scr . The parameters are fi = 3 and Q = 0.1 and e — 0.3 
(top row), e = 0.64 (bottom left), and e = 0.7 (bottom right). 



state. However, in contrast to the ID case where this colli- 
sion happens with the first excited state, in 2D the collision 
occurs with the second excited state, due to the absence of net 
topological charge in such a vortex-dipole (see top-right red 
curve) bearing two opposite charge vortices emerging from 
the central dip of the TF background. At this critical point 
e = Scr the dipole branch is unstable, having been destabi- 
lized through an SSB bifurcation at an e = e^cr > value 
(below which for e > the dipole is stable — see portion 
of red solid line in the figure). As this branch is followed 
(from top to bottom in the figure) a series of bifurcations oc- 
cur where the existing vortices are drawn to the periphery of 
the cloud, a dip in the center deepens leading eventually to 
a new vortex pair emerging (i.e., a higher excited state). In 
this manner the branches with even number of vortices are all 
connected. As more and more vortex pairs emerge, the cloud 
"saturates" and can no longer fit in new vortex pairs finally 
colliding with a dark soliton stripe (see lower blue branch in 
the figure). This overall bifurcating structure of even vortex 
numbers — with a s — >• — e symmetry where the solutions are 
just flipped by (x, y) — >> (— x, —y) — is depicted, with density 
and phase profiles, in the series of panels of Fig.HJb). 

As for the bifurcation scenario of odd number of vortices, 
the first excited state bearing a single vortex at the origin (for 
e = 0) is stable for small values of e, while it again sustains an 
SSB bifurcation for larger e. As e increases the vortex moves 
towards the periphery of the cloud and a dip at the center of 
the cloud deepens until a vortex pair emerges from it. This 
scenario connects the one- vortex branch with the asymmetric 

three-vortex (H h vortex tripole) branch, as it is depicted 

with the top (magenta and green) lines in panel HJa) and the 
series of snapshots in panels^c). As it is evident from the fig- 
ure, the asymmetric three-vortex branch eventually connects 
with the symmetric one for values of e —> 0. A similar bifur- 
cation occurs with the symmetric three- vortex branch, which 
becomes asymmetric with a deepening dip at the center where 
a vortex pair emerges (at the same time that a vortex is lost 
at the periphery), connecting in this way with the four- vortex 
branch [see series of snapshots in panels Htd)]. 
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FIG. 4: (Color online) The 2D generalization, (a) Bifurcation dia- 
gram for the 2D stationary nonlinear (vortex and DS stripe) states. 
Stable (unstable) branches, as per the corresponding BdG analysis, 
are depicted with solid (dashed) lines, (b) Series of density (left) and 
phase (right) configurations along the branch with even number of 
vortices corresponding to the circles in panel (a) [from top to bot- 
tom], (c) Same as (b) for the branch starting with one vortex and 
connecting with three vortices corresponding to the squares in panel 
(a) [from top to bottom], (d) Same as (b) for the branch starting with 
three symmetric vortices and ending with four vortices correspond- 
ing to the triangles in panel (a) [from top to bottom]. Parameter val- 
ues are: \i — 2 and Q — 0.2. The field of view for the configurations 
is [-10.5,10.5] x [-10.5,10.5]. 



As for the dynamics of unstable steady states, we have ob- 
served — in analogy with the ID case — that (a) a single vortex 
tends to migrate towards the minimum of the lossy side of the 
potential, while the remaining vortices (if present) perform al- 
most circular orbits at the periphery of the cloud where they 

(2) 

are eventually absorbed; and (b) past e = Sc r , using as an ini- 
tial condition the form of the TF background, also produces 



the spontaneous formation of an increasing number of vor- 
tices for larger values of e (namely, a "vortex sprinkler"). It is 
worth mentioning that the precise structure of the bifurcation 
diagram depends of the values of the propagation constant /i 
and the trap strength Q. For weaker ft and/or larger /i the ex- 
tent of the TF background will be larger allowing for a longer 
bifurcating chain of higher-order vortex states. Nonetheless, 
the displayed SSB instabilities and phenomenology and the 
nonlinear VT transition involving the cascade of blue-sky bi- 
furcations (notice that in the 2D case the order is reversed and 
the largest e bifurcation is that involving the TF and the dipole 
states) appear to be universal in confining PT-symmetric po- 
tentials. 

Conclusions. In the present work, we have developed some 
fundamental insights stemming from the interplay of defocus- 
ing nonlinearity and T^T-symmetric confining potentials. We 
identified both a symmetry-breaking bifurcation destabilizing 
the dark solitons that leads to non- stationary dynamics, as well 
as a nonlinear analogue of the VT transition that eventually 
terminates both the ground state and the dark soliton branch, 
yielding purely gain-loss dynamics within the system. Similar 
bifurcation phenomena and dynamics of mobility or of spon- 
taneous emergence of dynamical patterns forming out of the 
destabilization of the nonlinear states were identified in two- 
dimensional settings, for vortices. These investigations, we 
believe, pave the way for studying PT-symmetric systems 
in the context of defocusing nonlinearities and of higher di- 
mensional systems, which are some of the natural extensions 
of the 'PT-symmetric literature. A canonical set of investi- 
gations which is still missing concerns the effects of such po- 
tentials in three-dimensional continuum or higher dimensional 
lattice contexts, as well as the manipulation of nonlinear states 
emerging in these systems. These will be pursued in future 
works. 
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